On the inner Double-Resonance Raman scattering process in bilayer graphene 
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The dispersion of phonons and the electronic structure of graphene systems can be obtained 
experimentally from the double-resonance (DR) Raman features by varying the excitation laser 
energy. In a previous resonance Raman investigation of graphene, the electronic structure was 
analyzed in the framework of the Slonczewski-Weiss-McClure (SWM) model, considering the outer 
DR process. In this work we analyze the data considering the inner DR process, and obtain SWM 
parameters that are in better agreement with those obtained from other experimental techniques. 
This result possibly shows that there is still a fundamental open question concerning the double 
resonance process in graphene systems. 

PACS numbers: 63.20.D-, 63.20.kd, 63. 22. Re, 73.22.Pr 



In recent years, the physics of monolayer graphene has 
been thoroughly investigated, unveiling a wealth of in- 
teresting and unusual properties, most of which are re- 
lated to graphene's distinct electronic properties, that 
consist of a linear and isotropic dispersion of the elec- 
tronic states around the Fermi level (Ep) near the K 
point in the Brillouin zone (BZ). Bilayer graphene is also 
a very interesting material. While in the unbiased bilayer 
the valence and conduction bands touch each other at the 
Fermi level, a gap can be opened and tuned, for exam- 
ple, by the application of an external electric field [lJ-Q , 
which makes this a promising system for the fabrication 
of nanoelectronic devices. The development of bilayer- 
graphene-based bulk devices depends on the detailed un- 
derstanding of its electronic properties. Since the unit 
cell of AB stacked bilayer graphene is the same as that 
of graphite, one can model the bilayer electronic struc- 
ture using a tight-binding (TB) model for graphite [7|, by 
adapting the Slonczewski-Weiss-McClure (SWM) param- 
eterization la, |9fl of relevant couplings. There are several 
theoretical |10l - [l2| and experimental 0, [l3l - [T7| studies 
of theses SWM parameters, but the agreement between 
the reported values, obtained with different experimental 
techniques, is not entirely satisfactory. 

In previous resonance Raman studies of bilayer 
graphene [HI, [lj| in our group, the dispersion of the 
G' Raman band (also called 2D band) as a function of 
the laser energy was measured, and the nearest-neighbor 
hopping parameters 70, 71, 73 and 74 (shown in Fig. [T]) 




FIG. 1: The intra- (70 and t') and inter-layer (71, 73 and 74) 
tight-binding parameters in bilayer graphene. 



were determined. In Refllq. the fitting included also the 
in-plane second-neighbor hopping parameter t', which is 
expected to be of the same order as the out-of-plane 
nearest-neighbor parameters. The parameter A, which 
represents the difference between the on-site energies of 
the sublattices A and B, was also taken into account. 

Group theory analysis for bilayer graphene predicts 
four distinct DR processes (Pn, P22, P12, and P21) along 
the TKM direction, which are illustrated in Fig. [2] The 
triangularly-shaped isoelectronic curves around the K 
and K' points in Figs. [5Jc) and [5Jf ) are the equienergy 
contours of the 7r electrons involved in the scattering pro- 
cess. The value of the equienergy determined by the laser 
energy El, that creates the electron- hole excitations is 
the first step of the DR process, as shown in Figs. EUa-b) 
and[2jd-e). The analysis of the data in Refs. [l5| and Ha 
was done considering a double resonance (DR) process 
involving only backscattering of electronic states along 
the KT line with phonons along the KM line, a process 
we call the outer DR Raman process. From this analy- 
sis, we were able to obtain the SWM parameters, but our 
values for the 71 and 73 were found to be smaller than 
those determined using other experimental techniques. 

This restricted one-dimensional analysis of the DR pro- 
cess rests on the following assumptions: (i) that a one- 
dimensional integration along the TKM direction cap- 
tures the essential features of the DR process, as found 
in Ref. [HI; (ii) that some of the graphically determined 
double-resonant q vectors, related to forward-scattering 
processes [connecting points on the electronic equienergy 
curves surrounding points K and K' in Fig. Hfc)1, vanish 
by destructive interference, as found also in Ref. [lg (iii) 
that, by plotting the phonon density of states (PDOS) of 
graphene satisfying the DR process [19( , one can identify 
a strong singularity at the phonon q-vector involved in 
the outer process, and a much smaller PDOS value for 
the phonon q-vector along the KT line that backscatters 
an electronic state along the KM line, and also meets 
the DR condition, in a process we refer to as the in- 
ner process. This is the reason why attention is us ually 
paid solely to the one-dimensional outer process |19l - [21j . 
While the analysis in Ref. [ID did consider both outer and 
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inner processes, and the calculation of the Raman cross 
section considered all possible resonant and non-resonant 
processes, a critical approximation was employed: the 
matrix elements involved, related to electron-photon and 
electron-phonon couplings, were assumed to be constant 
and independent of the wavevectors k and q of the elec- 
trons and phonons, respectively. 

This state of affairs indicates that some of the con- 
ventional wisdom related to the DR process in graphene- 
based systems needs to be reevaluated. In the present 
work, we address the particular issue of the preponder- 
ance of the outer process. This may be particularly im- 
portant in the context of symmetry-breaking potentials 
imposed on the graphene system. For example, if a com- 
pressive or tensile strain is applied to the system, cither 
intentionally or due to interaction with a substrate, one 
expects the doubly degenerate G' peaks to split into two 
sub-bands G' + and G'_ , due to the movement of the Dirac 
cones and consequently the breaking of the symmetry of 
the DR Raman process. We reanalyze the data reported 
in Ref. [Tfj] considering now the inner DR process, in which 
the electronic excitations involved are along the KM line 
and the phonons are along the KT line. As discussed be- 
low, we obtain a very good fit of the experimental results 
using the same values for the parameters 70, 74, f, and A 
as in the outer-process case, but with values for 71 and 73 
which are in better agreement with those obtained from 
other experimental techniques. 

In order to understand the experimental results we 
obtained, it is important to build a bridge between ex- 
periment and theory. To achieve this, we must find a 
relation between the electronic and the phonon disper- 
sions of bilayer graphene. The electronic dispersion of 
bilayer graphene can be described in terms of the stan- 
dard SWM model for graphite, as mentioned above, using 
a TB model 0- Along the KT direction, this amounts 
to replacing the parameter a in McClure's expressions [1] 
by: cr = 70 [2 cos (2tt/3 - V3ka/2) + l] . Here k is the 
modulus of the electronic k- vector measured from the K 
point and a = 1.42 A is the in plane nearest- neighbor 
carbon distance. The bands in the bilayer are obtained 
from a TB Hamiltonian using the parameters 70 , 71 , 73 , 
74, t', and A. Along the high symmetry Kr direction, 
the 4x4 matrix factorizes and the dispersion of the four 
bands are given by: 



R 1 — 



-71 - <JV 3 + A' ± £1 



R 2 — 



71 + crv 3 + A' ± £ 2 



(1) 



where v, = 7^/70 , 
A' = A+t' 2 cos (2vr/3 - VZka/2^ +cos Ux/3- V3ka 
and 

£1 = yjili + vz<J ± A') 2 + 4((1 T vl)a* T ctv 3 (A' ± 7l )). 



It should be noted that these expressions are similar to 
those found in Ref. [15J , except that here we take A into 
account, to allow for the possibility of different doping 
levels between the two layers, and also include a t' pa- 
rameter for second-nearest-neighbor interactions within 
the same layer. For any of the Py (i, j = 1, 2) processes, 
we seek the dependence of the phonon energy E p h with 
El. In the initial step of this process (electron- hole cre- 
ation), the incident photon is in resonance with the ex- 
citation of the electronic state from the valence to the 
conduction bands at the kj point. In the following, we 
drop the vectorial notation for the k- and q- vectors, since 
we are considering only the TKM direction. The laser 
energy can then be written as: 



Er 



E^{h)-E ni {ki) , 



(2) 



which allows us to determine the momentum ki of the ex- 
cited electron in the process. The electron is then scat- 
tered from a state in the vicinity of the K point to a 
state in the vicinity of the K' point by emitting an iTO 
phonon with energy 

(3) 



E^M=E < {k i )-E^{k' j ) 



where qij depends on ki and kj. This equation uniquely 
determines the momentum k'j of the scattered electron, 

provided that E l A(qij) is known. The phonon energy 
can be computed, and is directly related to the Raman 
shift for a specific process, obtained with a given El. 
Physically, the difference between the outer and inner 
processes lies in the phonon wavevector qtj . 

As can be inferred from the geometry in Figs. [5Jc) 
and[2jf), qij — ki + kj for the outer process, and q^ = 
— (ki + kj) for the inner process, measured from the K 
point. In both cases, q^ has a maximum amplitude of 
approximately 2fc^. The outer and inner processes have 
their vectors ki and kj pointing in opposite directions. 
Interestingly enough, this means that calculations for the 
inner process can be done simply by switching ki into —ki 
and k'j into — kj in Eq. (1) and looking for values of kj 
satisfying Eq. (3). 

Figure [3] shows the experimental data (dots) and the 
TB fitting results (curves) for both the outer and the 
inner processes. Figure [3ja)outer shows the fit con- 
sidering the outer process for the DR scattering, i.e., 
with phonons along the KM line. The values obtained 
for 70, 71, 73, 74, A, and t' are shown in Table Ufa). 
Note that our value of 71 = 0.35 eV is slightly smaller 
than that usually found in the literature 71 « 0.40 
eV [HI, EH, HH ) from other techniques and calculations. 
However, the major discrepancy is in the value of 73 = 0.1 
eV compared to 0.30 eV found in the graphite litera- 
ture O HI, HI)]. Recent infrared studies in exfoliated 
bilayer graphene consider 73 = 0.30 eV [14J, but this 
value is not extracted directly from the experiments. 

In order to theoretically compute results and compare 
with the experimental data in Fig. 2, we have to con- 
sider a form for the phonon energy E l A(qij). In this pa- 
per, we consider a nonlinear relation for the iTO phonon 
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FIG. 2: (Color online) (a) Pn and P22 DR Raman outer processes involving the symmetric phonon. (b) P12 and P21 DR 
Raman outer processes involving the anti-symmetric phonon. (c) Outer process involving electrons on the Kr line and phonons 
on the KM line, (d) Pn and P22 DR Raman inner processes involving the symmetric phonon. (e) P12 and P21 DR Raman 
inner processes involving the anti-symmetric phonon. (f) Inner process involving electrons on the KM line and phonons on the 
Kr line. 
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FIG. 3: Laser energy dependence of the G' band peaks, with 
four different fittings, considering: ((a) outer) the outer scat- 
tering process with the best SWM parameters for this case, in 
Table |Ja); ((b) outer) the outer process with SWM param- 
eters of Table |T|b); ((b) inner) the inner scattering process 
with the best SWM parameters for this case, shown in Table 
|IJb); ((a) inner) the inner process with SWM parameters of 
Table |Ja) . Experimental data taken from Refflq. 



dispersion given by a second-order polynomial uu(q) = 
A + Bq + Gq 2 with q = ki + k'j and q = -(ki + fcj) for 
the outer and inner processes respectively, as explained 
above. We also use two distinct non-linear phonon dis- 
persions for the symmetric (S) and anti-symmetric (AS) 
branches. The best fit was obtained when we consider 
different dispersions for the two iTO phonon branches 
of bilayer graphene. Table [TT] shows the parameters ob- 
tained for the S and AS phonon branches, for both scat- 
tering processes. Note that the S phonon branch is more 
sensitive to the change of scattering process than the AS 
mode. 

Note also that, in Fig. [31 the fittings considering the 



TABLE I: Best values of the SWM parameters (in units of 
eV) obtained for (a) the outer and (b) the inner scattering 
processes. 





7o 


7i 


73 


74 


A 


t' 


(a) 


3.0 


0.35 


0.10 


0.18 


0.03 


0.10 


(b) 


3.0 


0.40 


0.30 


0.18 


0.03 


0.10 



TABLE II: Best values of the iTO phonon dispersion param- 
eters (w(q) — A + Bq + Cq 2 ) obtained for (a) the outer and 
(b) the inner scattering processes. 





symmetric 


anti-symmetric 




ABC 


ABC 




(meV) (meVA) (meVA 2 ) 


(meV) (meVA) (meVA 2 ) 


(a) 


146.3 85.6 -70.3 


150.4 65.4 -44.8 


(b) 


150 -65.9 -45.6 


153.4 -54.9 -37 



two different DR scattering processes produce rather dif- 
ferent results. In the inner-process case we cannot fit the 
data with the same parameters that produce the best 
fit for the outer process, shown in Fig. [3Ja)outer. This 
is shown in Fig. [3Ja)inner, where the theoretical curve 
obtained using the values of 71 and 73 from Table 1(a) 
differs considerably from the experimental values. The 
best fit for the inner process is shown in Fig. |3Jb)inner, 
with the values of 71 = 0.40 eV and 73 = 0.30 eV found 
in the literature, and with the same values of 7, A and t' 
as for the outer process. In a similar manner, the outer 
process cannot be fitted using the best-fit parameters for 
the inner-process case: deviations from the experimental 
results are again observed, as shown in Fig. [3jb)outer. 
The above discussion shows that the inner DR pro- 
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cess gives TB parameters in better agreement with those 
obtained by other experimental techniques, despite the 
fact that the number of phonon states with wavevectors 
that satisfy the DR process with visible photons is much 
larger for the outer DR process [18J. This result is pos- 
sibly an evidence that there is still a fundamental open 
question concerning the double resonance selection rules 
in graphene. 

In principle, there is a large number of phonon 
wavevectors that satisfy the DR process, connecting elec- 
tronic states in equienergy contours around the K and K' 
points. The calculation of the DR Raman profile must 
take into account all these possible DR phonon wavevec- 
tors, weighted by the q- vector dependence of the asso- 
ciated electron-phonon matrix elements. Some interfer- 
ence effects also appear when the Raman expression is 
squared out in order to calculate the Raman intensity 
[241 ]. To our knowledge, there is only one full calcula- 
tion of the shape of the Raman DR bands that takes into 
account the q-vector dependence of the electron-phonon 
matrix elements, and this calculation predicts an asym- 
metric shape for the DR bands in monolayer graphene 
[24j . However, the experimental DR Raman bands ob- 
tained with visible photons bands are nicely fitted by a 
single Lorentzian lines. This possibly shows that some 
ingredients are missing in order to fully explain the DR 
process in graphene systems. 

The results obtained in this paper could be the starting 
point to investigate other systems which constitute a hot 
subject in graphene physics, such as strained or twisted 
bilayer graphene. In the former case, bilayer graphene 
is grown on an insulating material (such as Si02), which 



imposes a strain on the graphene system. This setup has 
been considered as a building block for microelectronics 
[25l ]. In the latter case, the stacking of the two layers 
is different from the usual AB (Bcrnal) stacking, being 
similar to what is found in naturally occurring and syn- 
thetic crystals presenting a variety of defects, which affect 
the stacking order mainly in the c-axis direction [26l . l27j . 
In both cases, there is a significant modification in the 
electronic and optical properties of the systems, which 
directly influences the Double Raman bands. This pro- 
duces different G' peaks that can be theoretically studied 
by looking at the P^- processes. 

Finally we hope that future experimental work on 
graphene systems will reveal more insights on the in- 
ner process and on scenarios where it can be considered 
more dominant than its counterpart, the outer process. It 
would also be quite fascinating if one could theoretically 
calculate the Raman cross section of bilayer graphene 
or of any of its exotic descendants such as twisted or 
strained graphene systems by integrating over all k points 
on the isoenergetic surfaces without any approximations. 
This might show that points other than those along the 
KT and KM lines need to be taken into account, which 
might lead to a better understanding of the inner and 
outer processes. 
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em Nanotubos de Carbono - MCT, and the Brazilian 
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ter for Integrated Nanotechnologies, jointly operated by 
Los Alamos and Sandia National Laboratories. 
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